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I. INTRODUCTION 



II. CONTRIBUTION OF THE HYPERFINE 
INTERACTIONS 



The quark model predicts that the masses of the 
mesons and baryons are given by the sum of the con- 
stituent masses of the quarks and the hyperfine splitting 
[IHIj]. As described in [l|, the term with hyperfine split- 
ting emanates from a spin-spin interaction produced by 
one gluon exchange, which also depends on the masses of 
the quarks (for mesons this term contains a mass sum of 
a quark pair and for baryons it contains a sum over all 
the possible quark pairs). 

In Refs. [H the baryons were studied in a quark model 
with color hyperfine interactions. Some parameters of 
meson and baryon spectroscopy, including some predic- 
tions of the quark masses, can be found in In Ref. 
Q, from the meson and baryon relations, the authors 
fitted both the mass differences and mass ratios with a 
single set of quark masses. In Q estimations of masses 
were made via magnetic moments of baryons and via the 
vector meson masses. In Refs. [8, 9] it was shown that 
in the constituent quark model, the Feynman-Hellmann 
theorem and semiempirical mass formulas can be applied 
to give useful information about the masses of mesons 
and baryons. 

In this paper we use a schematic study (two quark 
interaction) to calculate meson and baryon theoretical 
masses from which the constituent quark masses can be 
obtained. This paper is organized in the following way: 
in Section II we give complete method how to calculate 
theoretical masses, in Section III we explain fitting pro- 
cedure and our main results, discussion and conclusions 
are given in Section IV. 



Using the colored version of the Fermi-Breit hyperfine 
interaction (FB HFI) [l0rlT3| and the Glozman-Riska hy- 
perfine interaction (GR HFI) [IM!], we present estima- 
tions of the theoretical meson and baryon masses. Strong 
FB HFI Hamiltonian [l3| has the following form: 



Hfb — C 



i>j 



(1) 



and it has explicit color and spin exchange dependence 
and implicit (by way of quark masses) flavor dependence. 
Here &i are the Pauli spin matrices, Af are the color Gell- 
Mann matrices and C is a constant. FB contribution to 
hadron masses is given by: m„FB = {v\ (xl ^fb |x) \ v )i 
where \ denotes the spin wave function and v - the flavor 
wave function. Strong GR Hamiltonian [14| is of the 
form: 



#GR = -G x £ ("IP (Af Xf ) (; 



(2) 
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where Af are Gell-Mann matrices for flavor SU(3), are 
the Pauli spin matrices and C x is a constant. We employ 
this schematic flavor-spin interaction between quarks and 
antiquarks which leads to Glozman-Riska HFI contri- 
bution to hadron masses: to^gr = (v\ (x| -Hgr |x) \ v )i 
where m, are the constituent quark effective masses: m u 
— fid 7^ m s and v - flavor wave functions. 

For mesons, there are two flavor SU(3) multiplets ac- 
cording to product: 3 ® 3 = 1 + 8, i.e. one singlet and 
one octet. For baryons, according to: 3 <E) 3 ® 3 = 
1 + 8ms + 8 ma + 10, there is one singlet, two octets 
(one mixed symmetric (MS) and one mixed antisymmet- 
ric (MA)) and one decuplet. Young diagrams for these 
SU(3)f multiplets, as well as the weight diagrams, can be 
found in textbooks (see e.g. Regarding baryons, 

in case of octet members (as they have the MS and MA 
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part of the flavor wave function), the contribution of 
HFI is calculated as sum of the symmetric and antisym- 
metric part, so we start from the following expression: 
\v t) = 71 Wms + v M a)- 



III. FITTING PROCEDURE 

For our calculations we used multidimensional least- 
square fit of masses, using subroutine "lfit" from Numer- 
ical Recipes in FORTRAN 19], modified according to 
the instructions in the last paragraph of §15.4 of [19| . 
We fitted simultaneously all parameters: m u , m s , m c 
and the constant C in such a way that we minimized x 2 
between the measured and theoretical masses. Equations 
for meson and baryon theoretical masses (see §IV) are at 
first linearized by expansion in Taylor's series up to the 
first order, and in that way we obtained the correspond- 
ing system of linear equations for differences between ex- 
perimental and theoretical masses in which the unknown 
variables are corrections of the parameters. These correc- 
tions, as well as their corresponding uncertainties, are de- 
termined using linear least square method, and then the 
parameters are accounted for the values of these obtained 
corrections. With these new (corrected) parameters the 
previous procedure is repeated until the fit converges, 
that is while the \ 2 value between the experimental and 
theoretical masses decreases. In that way, after several 
iterations, we obtain final values for the parameters and 
their uncertainties. The uncertainties are estimated dur- 
ing fitting procedure as square roots of the corresponding 
diagonal elements of covariance matrix, according to eq. 
(15.4.15) of Ref. [H. 

For every analyzed system of equations a fast conver- 
gence is obtained, even in the case when starting values 
of parameters differ very much from their final values, 
which favors the correctness of our theoretical model as 
well as our fitting method. 



IV. RESULTS AND DISCUSSION 

We calculated FB and GR contributions to meson 
and baryon masses. The mass formulas of mesons and 
baryons with FB HFI included are given in eqs. ((3]) - (j4}, 
while their mass formulas with GR HFI included can be 
found in [l5| (eqs. (6)-(10) for mesons and eqs. (11)-(13) 
for baryons). 

Here we study the hadron masses using FB and GR 
HFI in schematic approximation (two-particle interac- 
tion). The masses of constituent quarks m u (= m^), 
m s , m c and the constants C m and C b are calculated from 
least-square fit of the theoretical equations for meson and 
baryon masses with HFI included. The corresponding ex- 
perimental masses are taken from " Particle Data Group" 
site: |http://pdg.lbl.gov| [20j . 

We give the equations for meson masses with FB in- 
teraction included for the following mesons: light pseu- 



doscalar mesons n, K. light vector mesons p, K*, u, <p, 
charmed mesons D, D* , strange charmed mesons D$, 
Dg and double charmed mesons r) c , 3 /if) (eqs. ©)• In 
these equations, the constant for mesons is denoted by 
C m . We did not calculate r\ and rf contribution because 
of their mixing and they cannot be described within such 
a model. Because of that mixing (due to the same quan- 
tum numbers) , their flavor wave functions are given only 
in a first approximation (2lj and therefore calculations 
are not sufficiently precise. The mixing of the states also 
changes the properties and shifts masses from the theo- 
retical predictions. 

3C m 3C m 
ni n = 2m u — ; itik = m u + m s 



mf. m u m s 

m„ = 2m u H 5- ; mg, = m u + m s H 

mi. m u m s 

C C m 

m u = 2m u H T ; m v = 2m s H T 

mf. mi 

3C m C r ' 
m D = m u + m c ; m Dsr = m u + m c + 



m„m r 

m Ds = m s + m c ; m Ds , = m s + m c + 



m„m r 



m nc = 2m c 



3C" 



m s m c 



2m r 



m s m c 



(3) 

We also present the theoretical mass equations with 
FB HFI included for the following baryons: light baryon 
octet N, S, S, A, light baryon decuplet A, £*, 3*, and 
heavy baryons S c , 3 CC , A c , S*, Q c (eqs. (HJ)). Constant 
for baryons is denoted by C . 



m N 



= 3m u - 3C fc 



1 



2m?, 



1 ( 1 



= 2m u + m s + 2C b — - 



1 / 1 



m s = m u + 2m s + 2C"^r - 



m . 



m , 



771A 

ma 

TOHcc 

m Ac ~- 
mn c - 



= 2m u + m s — 3C b 



2m?, 



= 3m M + 3C fc 



1 



2m? 



2m„ + m s + C b 



1 ( 1 



1 fl 



= m u + 2m s + C— - + — 



3m s + 3C* 



1 



rn f 



rn u 



(4) 



2m? 



= 2m u + m c + 2C l 



1 fl 



rn u + 2m c + 2C l 



2m u + m c — 3C b 



1 fl 



m£ \ 4 



rn , 



2m? 



1 fl 



= 2m u + m c + C b — - + — 



2m s + m c + C b -^ (\ 
mi \ 2 



m 
m 
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TABLE I: Constituent quark masses obtained from the meson fits. For each combination of mesons, the results are obtained 
from their masses with FB HFI included (upper rows) and with GR HFI (lower rows). Note that constant C m differs for the 
two HFIs. 



Fit 


Mesons 




Quark masses (MeV) 




r 


x 2 


No. 




m u = m d 


m 3 


m c 


(xl0 7 MeV 3 ) 


(xlO 6 ) 


1 


TT, K, p, K* , (J, If, D, 


318.12 ± 0.11 


477.29 ± 0.21 


1591.42 ± 0.21 


1.5977 ± 0.0017 


0.095336 




£>*, D S , D* s , 77c, J/</> 


304.55 ± 0.11 


534.33 ± 0.19 


1576.94 ± 0.21 


2.1331 ± 0.2346 


0.304832 


2 


tt, K, p, K* , to, p 


308.05 ± 0.12 


484.91 ± 0.24 




1.5043 ± 0.0017 


0.000411 






292.35 ± 0.12 


552.80 ± 0.22 




2.0284 ± 0.0023 


0.143289 


3 


D, D* , D S , D* s , 


454.75 ± 0.47 


547.48 ± 0.49 


1524.36 ± 0.35 


2.6936 ± 0.0095 


0.003475 






453.68 ± 0.44 


546.82 ± 0.49 


1524.94 ± 0.35 


4.0280 ± 0.0138 


0.003480 



TABLE II: The same as Table I, but for baryon fits. Note that for heavy baryons fit did not converge and constitutive quark 
masses could not be obtained. 



Fit 




Baryons 




Quark masses (MeV) 




C b 


x 2 


No. 






m u = m d 


m s 


m c 


(xl0 7 MeV 3 ) 


(xlO 6 ) 


1 


N, E, 


S, A, A, E*, H*, 


381.53 ± 0.09 


537.26 ± 0.18 


1368.15 ± 0.27 


1.3454 ± 0.0021 


0.153491 




n, E c 


j ^cc; Ac, S c , r^c 


541.04 ± 0.18 


655.00 ± 0.22 


1358.62 ± 0.25 


2.4878 ± 0.0047 


0.125597 


2 


N, E, 


S, A, A, E*, H*, 


362.52 ± 0.10 


538.68 ± 0.20 




1.2789 ± 0.0020 


0.001165 








500.42 ± 0.20 


621.20 ± 0.25 




1.6951 ± 0.0043 


0.002719 



TABLE III: Constituent quark masses obtained from hadron fits when hadrons with two c-quarks are excluded. The upper 
rows correspond to FB HFI, and the lower rows to GR HFI. 





Hadrons 


m u = m d 


Quark masses (MeV) 
m s 


m c 


C 

(xl0 7 MeV 3 ) 


x 2 

(xlO 6 ) 


Mesons 


tt, K, p, K* , u, p, D, 


314.75 ± 0.12 


466.80 ± 0.21 


1627.31 ± 0.25 


1.5546 ± 0.0017 


0.033123 




D*, D s , D* s , A c , E*, Q c 


300.88 ± 0.11 


526.65 ± 0.19 


1605.62 ± 0.26 


2.0660 ± 0.0023 


0.263221 


Baryons 


N, E, S, A, A, E*, H *, 


365.69 ± 0.09 


530.08 ± 0.18 


1700.17 ± 0.37 


1.2513 ± 0.0019 


0.015451 




E c , A c , E c , f2 c 


506.26 ± 0.19 


623.40 ± 0.22 


1629.69 ± 0.35 


1.8163 ± 0.0041 


0.032378 



TABLE IV: Predicted masses (in MeV) of constituent quarks in different phenomenological models. 



Mass 


Ref. [5] 


Ref. [8J 


Ref. [1] 


Ref. [22J 


Ref. [6J 


Ref. \7} 


Ref. [15] 


m u = m d 


220 


300 


310 


290 


360 


337.5 


311 


m s 


419 


475 


483 


460 


540 


486 


487 


m c 


1628 


1640 




1650 


1710 


1550 


1592 



For each set of equations, the minimized x 2 values 
for masses are calculated by the formula: \ 2 — 
12f=i ; where Tj is the model prediction for the 

hadron mass, Ei the experimental hadron mass and at 
the uncertainty of the mass. 

The values of constituent quark masses and their un- 
certainties (in MeV) are given in Tables HI and HI1 It is no- 
ticeable that the fits are satisfactory for both HFIs and in 
cases of all mesons, light mesons and light baryons. They 



show that the constituent mass is only slightly modified 
by the dynamics of confinement. The constituent masses 
differ for the quarks in mesons and those in baryons. This 
is not unexpected since they are confined in quite differ- 
ent systems. Also, as noted in [§| there are no theoreti- 
cal reasons why the masses determined from the baryon 
fits should coincide exactly with those determined from 
the meson fits because these quark masses are constitu- 
tive masses, i.e. effective ones. Roncaglia et al. [1, 
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also showed that constraints (inequalities) for mass dif- 
ferences are stronger for baryons than for mesons, and 
therefore requirement for the same set of quark masses 
for both baryons and mesons would result in significantly 
poorer fits [8|. 

In the case of heavy mesons, the fit resulted in some- 
what increased values for u and s quarks, and in the case 
of heavy baryons fit even did not converge. Therefore, 
our results suggest that both of the analyzed HFIs can 
be used in a first approximation for modeling the term 
with hyperfine splitting in the quark model only for light 
mesons and light baryons. That is both HFIs well de- 
scribe systems which contain quarks from SU(3)f, i.e. 
light quarks: u, d, s. In order to verify this conclusion, 
we also show the results for the case when hadrons with 
two c-quarks are excluded (see Table lITij) : meson fit with- 
out r] c and J/tp-, and baryon fit without S cc . From these 
two fits we can calculate mass differences m s — m u and 
m c — m s which may be then compared with the corre- 
sponding inequalities given by Roncaglia et al. [1, Q, 
obtained from the Feynman-Hellmann theorem. We can 
see from Table HiTl that in this case the constituent quark 
masses obtained for FB HFI satisfy all of the inequalities 
given in dlH, contrary to Tables U and HI where FB HFI 
breaks some of these inequalities. Also, from Table IIIII 
one can see that GR HFI violates one or more of the in- 
equalities, but to a lesser extent than in Tables HI and [TT1 
Even more, in this case the obtained constituent quark 
masses are more realistic than those from Tables U and 
HTl From each set of equations (except for heavy baryons) 
we calculated masses of constituent quarks and constants 
C m and C b , and we got slightly different values, but all 
of them are from expected mass ranges of constituent 
quarks. As one can see from Table IIV1 our predictions 
for constituent quark masses are in the range of masses 
obtained using different phenomenological models. When 
comparing these two HFIs, it is interesting to note that 



the values of obtained quark masses are similar although 
one interaction is color-spin, and the other one is flavor- 
spin. As it can be seen from Tables U and HH two different 
HFIs result with masses of constituent quarks which are 
more similar than those obtained by the same HFIs ap- 
plied to different groups of hadrons (e.g. light mesons, 
heavy mesons, light baryons, ...). For example, results 
for m u (= wid) and m s for FB HFI applied on light and 
heavy mesons are more different than results for m u (= 
rrid) and m s for FB and GR HFIs applied on light mesons. 

When comparing FB, the color-spin interaction, with 
GR HFI, which is flavor-spin interaction, it is interesting 
that we obtained similar results. We can conclude that 
both HFIs describe masses of light mesons and baryons 
very well. It shows that both HFIs well describe systems 
which contains quarks from SU(3)f- In case of heavy 
baryons the fits for both HFIs did not converge, which 
might indicate that quark model which includes these two 
HFIs is not satisfactory approximation for heavy baryons. 
It might also indicate that FB and GR HFI are not the 
complete effective two-quark interactions and therefore 
they cannot be successfully applied to quark systems out 
of SU(3) group. If we do not take into account hadrons 
with two c-quarks the FB HFI becomes a good approxi- 
mation, even for hadrons having one heavy c-quark. We 
also show that the constituent quark masses are very sen- 
sitive to the system in which they are confined, and their 
values differ less or more in different systems (for exam- 
ple heavy mesons and light mesons). We obtained that 
least-square fit is more accurate for FB than for GR in- 
teraction in case of light mesons, but in case of heavy 
mesons and light baryons accuracy of the fit is similar. 
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